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Summary 

A progressive damage methodology has been developed to predict damage growth and residual 
strength of fiber-reinforced composite structure with through penetrations such as a slit. The methodology 
consists of a damage-dependent constitutive relationship based on continuum damage mechanics. Damage is 
modeled using volume averaged strain-like quantities known as internal state variables and is represented in the 
equilibrium equations as damage induced force vectors instead of the usual degradation and modification of the 
global stiffness matrix. 


Introduction 

The characterization of damage in notched composite laminates has led to the development of 
progressive damage models to predict damage initiation, damage growth, and residual strength. Allen and 
Harris developed a damage-dependent constitutive model [1,2] which utilizes kinematic-based volume 
averaged damage variables to represent the effects of matrix cracking and fiber fracture. This model has a 
matrix crack growth law for fatigue as well as monotonic tension and the kinematic effects of delaminations are 
modeled empirically. An experimental verification of the accuracy of the model to predict stiffness loss in 
toughened composite systems was documented by Coats [3]. The progressive damage analysis algorithms of 
the Allen-Harris progressive damage methodology were implemented into a general purpose finite element 
code developed by NASA, the Computational Structural Mechanics Testbed (COMET) [4]. The methodology 
was developed [5] to model mode I and n matrix cracking and tensile fiber fracture in composite center-notch 
panels that were loaded monotonically in tension. The numerical details of accounting for damage and 
calculating internal state variables are presented in this paper. The paper will conclude with a comparison of 
analytical predictions with experimental results for AS4/938 center-notch tension panels. 

Progressive Damage Analysis 

The progressive damage analysis consists of a constitutive model where internal state variables (ISV) 
represent the average effects of local deformation due to the various modes of microcrack damage. The 
constitutive model for an individual ply within a finite element may be written as [ 1] 

<^ij L = Qijki {eid - «ki} L (l) 

where aij L are the locally averaged components of stress, Qjjki are the elastic moduli in ply coordinates, and 

£kl L are locally averaged components of strain. The internal state variables, cxy^, represent the local 
deformation effects of the various modes of ply-level damage. The damage dependent constitutive equations 
(1) can be transformed into laminate equations and substituted into the laminate resultant force and moment 
equations. This results in damage-dependent laminate resultant force and moment equations and can be written 
as [4,6] 
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{N} = [A]{£°}+[B]{K L } + {fd} 

{M} = [B]{£^} + [D]{K L } + {gd} 

where {N} and {M} are the resultant force and moment vectors, respectively; [A], [B], and [D] are the 
laminate extensional, coupling, and bending stiffness matrices, respectively; {e£} is the midplane strain vector; 

{k l } is the midplane curvature vector; and {f d } and {g d } are the damage resultant force and moment vectors, 

respectively [4,7]. The application of {fd} and { g d } to the undamaged material will produce midplane strain 
and curvature contributions equivalent to those resulting from the damage-induced compliance increase. The 
damage resultant forces and moments are written as: 
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where [Q] k is the transformed elastic modulus matrix for the kth ply in laminate coordinates. When equations 

(2) and (3) are substituted into the plate equilibrium equations, the result is a set of governing differential 
equations which can be integrated against variations in the displacement components to produce a weak form 
of the damage-dependent laminated-plate equilibrium equations [4,6]. By substituting displacement 
interpolation functions into the weak form of the plate equilibrium equations and following well known finite 
element procedures, the assembled equilibrium equations are obtained as: 

[K] [8] = [F A ] + [F D ] (6) 

where [K] is the global stiffness matrix of the undamaged structure, [8] is the global displacement vector, and 
[Fa] is the applied force vector. The damage induced force vector, {Fd}, is calculated from the element 
shape functions and the damage resultant forces and moments. Since the effects of damage are represented as 
damage-induced force vectors on the right-hand side of equation (6), the element stiffness matrix does not need 
to be recalculated as damage progresses as long as the nonlinearity in the load-deflection curve is not large. 

The model requires only the standard ply stiffness and strength properties determined from 
unidirectional laminates for characterizing the material and predicting residual strength. The progression of 
damage is predicted by an iterative and incremental procedure outlined in the flow chart shown in Figure 1. 
This entire progressive failure analysis scheme has been implemented into the finite element formulation in the 
NASA Computational Mechanics Testbed (COMET) [4] computer code. The first block of Figure 1 is a 
description of the finite element model. Block numbers are shown in the right hand comer of each box in the 
flowchart. Blocks 2 and 3 are processors that calculate the element stiffness matrices and assemble and factor 
the global stiffness matrix. The compliance changes due to damage are accounted for by combining the 
damage induced force vector with the applied force vector and solving for the global displacements in equation 
(6). This solution process occurs in blocks 4 and 5 and then the element stress resultants are computed in 
block 6. 

The ply-level strains and stresses are computed in block 7 and as long as the strains in a material 
element (local volume element or finite element) are less than the critical strains, e k ] crit , no damage exists and 
the internal state variables have zero values. When the strains reach their critical value, the element is damaged 

and this damage is represented by an internal state variable. The formulation for a k i was derived based on the 
assumption that when fiber fracture, mode I (opening mode) matrix cracking, or mode II (shear mode) matrix 
cracking occur in a ply within an element, the longitudinal, shear, and transverse stresses for that plv in that 
element are J 
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where S Cf , S£ r , and S c ^ are the ply longitudinal, transverse, and shear critical strengths, respectively; y, p, 

and vjr are scale factors that describe the load carrying capability of the material after the occurrence of tensile 
fiber fracture, mode I matrix cracking, and mode II matrix cracking, respectively. This is illustrated in Figure 
2 where it was assumed that for mode I matrix cracking and tensile fiber fracture, 90% of the failure load in the 
damaged ply of the damaged element would be redistributed to adjacent plies and elements. This was 

accomplished by setting y and P to 10%. The load redistribution for mode II matrix cracking was modeled 



Figure 1. Progressive Failure Analysis. 


after Iosipescu shear test data [5, 8]. 
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Figure 2. Load Carrying Capability After Ply Failure. 


The formulation for aki was obtained by substituting equations (7), (8), and (9) into the damage- 
dependent constitutive equations (1) and applying the following equation 

aid = ajf 1 + da k i (10) 


where aid is the updated internal state variable and a£ ld is the internal state variable of the previous damage 
state. The resulting equations are 
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where dan, da 22 , and dan are the incremental changes in the internal state variable for tensile fiber fracture, 
mode I matrix cracking, and mode II matrix cracking, respectively; £n, £ 22 , and y 12 are the longitudinal, 
transverse and shear strains, respectively; and Q is the ply-level elastic modulus in ply coordinates. The 
damage state is updated in block 8 using equation (10) and the ply-level stresses and strains are post- 




processed. Note that damaged elements are not removed. Instead, the damage state in the damaged elements 
increases linearly with strain. 

The equilibrated solution is obtained when the damage-induced force vector and the change in 
displacements become negligible. If the solution shows that equilibrium is established, then the next load 
increment is applied. If equilibrium is not established, then iteration is performed holding the loads at the 
current level. As the applied load approaches the catastrophic failure load, the number of iterations required to 
establish equilibrium increases. This is because the strains in the surrounding undamaged elements approach 
critical values. In the process of iteration, displacement increases lead to increases in strain levels in these 
elements and the strains may exceed the failure strains. Therefore, damage may grow from one iteration to the 
next. Eventually, as loads increase, a load level will be reached where equilibrium is unattainable. This value 
of load is defined to be the catastrophic failure load, or residual strength. 

Results 

Center-notch tension panels, as shown in Figure 3, were loaded monotonically to failure. Linear 
elastic fracture mechanics (LEFM) predictions are compared to the progressive damage model predictions and 
the experimental data in Figure 3 for four notch sizes and three panel widths. The experimental data is 
displayed with error bars denoting the maximum and minimum values of the measured strength data. The 
LEFM prediction, which is plotted as the solid line, was obtained by determining the fracture toughness of the 
laminate from the residual strength of the smallest specimen with a (2a/W) ratio equal to 1/4. The experimental 
values for the 22.9 cm notch are about 45% higher than the LEFM based predictions, whereas the progressive 
damage model correctly predicts this increase in the fracture strength. From these results it is obvious that the 
finite element model predictions are more accurate than LEFM for wide panels where fracture resistance effects 
are dominant. Mesh refinement issues concerning the center-notch configuration are discussed in detail in 
reference [5]. 



Figure 3. Comparison of LEFM and Model Prediction for the AS4/938 Center-Notch Tension Panels. 
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Concluding Remarks 

A progressive damage methodology has been developed to predict damage growth and residual 
strength of fiber-reinforced composite structure with through penetrations such as a slit. The methodology 
consists of a damage-dependent constitutive relationship based on continuum damage mechanics. Damage is 
modeled using volume averaged strain-like quantities known as internal state variables. The progressive 
damage analysis algorithms have been implemented into a general purpose finite element code developed by 
NASA, the Computational Structural Mechanics Testbed (COMET). Strain failure criteria are used for damage 
initiation and the internal state variables are used to numerically simulate damage growth. Damage is 
represented in the equilibrium equations as damage induced force vectors instead of the usual degradation and 
modification of the global stiffness matrix. Load is applied incrementally and a classical iterative method is 
used to establish equilibrium. Analytical predictions of residual strength were made for composite panels with 
central through slits and the computational results agree within ± 10% of the experimental results. 
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